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We propose a new scheme to realize holonomic quantum computation with rf-SQUID qubits in 
a microwave cavity. In this scheme associated with the non-Abelian holonomies, the single-qubit 
gates and a two-qubit control-Phase gate as well as a control-NOT gate can be easily constructed 
by tuning adiabatically the Rabi frequencies of classical microwave pulses coupled to the SQUIDs. 
The fidelity of these gates is estimated to be possibly higher than 90 % with the current technology. 
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Since the proposal of holonomic quantum computation 
[1] , research on quantum gates based on Abelian or Non- 
Abelian geometric phase shifts has attracted significant 
interests both experimentally and theoretically [2-10]. It 
is believed that these quantum gates could be inherently 
robust against some local perturbations since the Abelian 
or Non-Abelian geometric phases(holonomies) depend 
only on the geometry of the path executed. On the other 
hand, quantum information processing using Josephson- 
j unction systems coupled through a microwave cavity has 
been paid particular attentions recently [11-19]. Never- 
theless, how to realize the holonomic quantum compu- 
tation using superconducting quantum interference de- 
vices(SQUIDs) in a cavity has not been addressed. 

In this paper, we propose a novel scheme to achieve 
holonomic quantum computation using SQUIDs in a cav- 
ity. Based on the non-Abelian holonomies, two non- 
commutating single-qubit gates and a two-qubit control- 
Phase gate as well as a control-NOT gate are realized 
by tuning adiabatically the Rabi frequencies of classical 
microwave pulses coupled to the SQUIDs. The distinct 
advantages of the present scheme may be summarized as 
follows, (i) The energy spectrum of each SQUID qubit 
may be adjusted by changing the bias field; (ii) the strong 
coupling limit g 2 > > ("yn) may be easily realized, where g 
is the coupling coefficients between the SQUID qubit and 
the cavity field, n the life time of the photon in the cav- 
ity and 7 the life time of the excited state of the SQUID 
qubit; (iii) the decoherence caused by the external envi- 
ronment can be significantly suppressed; (iv) the fidelity 
of these gates may be higher than 90 % with the current 
technology. 

We consider an rf-SQUID (with junction capacitance 
C and loop inductance L) in an microwave cavity. The 
Hamiltonian of the rf-SQUID can be written as [15,20] 

Hs = 9l + (*-**) 2 _ EjCOS ( 27r ±) (1) 

2C* 2L V *o/ 

where Ej is the maximum Josephson coupling energy, 
$ x the external magnetic flux and <&o = h/2e the flux 
quantum. The conjugate variables of this system are the 
total charge Q and the magnetic flux $ which satisfy 



[<&, Q] = ih. It is well known that the Hamiltonian of Eq. 
(1) is quite similar to that of a particle moving in a double 
well potential. By changing the device parameters C, L 
and the control parameters Ej, <$> x , one can control the 
structure of energy levels in the SQUID. 

Let us address a (3+l)-type system with three lowest 
levels (|a ), \a\), \g)) and an excited level (|e)) in the 
SQUID (see Fig. 1). In the system, the \g) <-> |e) transi- 
tion with energy level difference u> eg is coupled to an one- 
mode cavity field with frequency lo c and the \a,i) «-* |e) 
transition with the energy level difference Lu e i (i = 0,1) 
is coupled to the classical microwave pulse with the mag- 
netic component as Bi (r, t) cos (toot) where LOi is the en- 
ergy difference between the states |a^) and |e) and can 
be adiabatically changed. We may ensure that the " 3- 
photon resonance" condition, i.e. uj eg — lu c — uj e i—uji = A, 
is satisfied. In the interaction picture, the Hamiltonian 
of the system can be written as 

H = A |e) (e| + ft (t) |e> (a \ + r>i (t) |e) (a x | (2) 
+ga\e) (g\ + h.c. 

where a (a^) is the photon annihilation (creation) opera- 
tor of the cavity field. Here the coupling coefficients g can 
be written as g — L^ 1 \ju c j (2/xoH) (g\ $ |e) f s B c (r) • dS 
where fiQ is the vacuum permeability, S the surface 
bounded by the SQUID ring and B c (r) the magnetic 
component of the cavity field [15]. The Rabi frequencies 
fli (t) (i = 0, 1) are proportional to the matrix elements 
(ai \ $ |e) [15]. It is pointed out that the detuning A of 
the SQUID can be adjusted by changing the bias field. 




1 



FIG. 1. A schematic diagram of the energy level in the 
SQUID coupled to the single mode cavity field (with cou- 
pling constant g) and two microwave pulses (with coupling 
constants Qo and ill). The 3-photon resonance condition is 
satisfied and A is the detuning. 

Regarding single-qubit operations, our scheme is sim- 
ilar to that proposed in [4,5]. We choose the states |oo) 
and \ai) as the computational basis and \g) as an an- 
cillary state. When Oo = ^1 = 0, the states |aj) |0) c 
(i = 0, 1) span an eigenspace of the Hamiltonian (2) 
with zero eigenvalue., where, |0) c is the vacuum state 
of the cavity field. If the cavity is cooled to zero temper- 
ature and the quantum gate operations are switched off, 
i.e., the Rabi frequencies of all the classical microwave 
pulses are set to zero, the state a\ao) + (3\ai) of the 
qubit is isolated from the state of the cavity and does 
not change with time. When the Rabi frequencies fio 
and Sli change adiabatically along a close path C in the 
parameter space M and return the point correspond- 
ing to fl = ill = (we refer to this point as O), an 
initial state |^o) = ( a \ a o) + P \ a i)) |0) c of the qubit- 
cavity composite system evolves according to the rule 
|*o) -> f/(C)|* ) [5,21]. Here, U (C) = Pcxp J c A is 
the non-Abelian holonomy associated with the path C 
and A = ^2 A^dX^ is the U (2)-valued connection ex- 
pressed as A % i = (Di (A)| gj- \Dj (A)) where {A M } are 
the coordinates of the parameter space M and \D t (A)) 
(i = 0, 1) the basis of the eigenspace of the Hamilto- 
nian (2) with zero eignevalue (hereafter we refer to the 
basis as dark states). Since the holonomy U (C) is a uni- 
tary transformation in the space spanned by the states 
| ao) |0) and |ai) |0) c , it can actually be considered as a 
unitary transformation that only acts on the qubit state 
which is the superposition of \ao) and \a\). 

Without loss of generality, we assume that the cou- 
pling coefficient g is real and positive, and choose Qo = 
3 tan (^e^ and Qi = g tan (9) sec (£) where 0,£ e 
[0, f) and <j>i e [0, 2?r) (i = 0, 1). We take the angles 
£, 9, (fro and tfri as the coordinates of the parameter 
space M. The dark states of this invariant subspace 
spanned by the states {|aj) |0) , \g) |1) , |e) |0) } of the 
Hamiltonian (2) can be written as the vector functions 
in M: \D ) = cos (f) \a Q ) |0) c - sin (f) e**° \g) |l) c and 
|£>i) = -sin(6»)sin(C)e 1 ^ 1 "^) |o ) |0) c - 

sin (0) cos (0 |g) |l) c + cos (0) |oi) |0) c . We wish to 
point out that our choices of the coordinates and dark 
states are quite different from those in Refs. [4,5]. Here, 
the dark states \D ) and \D\) are single-valued at the 
point O(0 = £ = 0), e.g. \D t (0)> - k) |0) c . 

It is well known that any single-qubit gate operation 
can be decomposed into the product of rotations about 
axies z and y: R z (ip) = e lipa * and R y (ip) = e 4 ^ , 
where, ip is the angle, a z , and o~ y are Pauli matri- 
ces defined as a y = i (|ai) (ao| — |ao) and ct z = 
(|ao) (ao| — |ai) (ai |). Therefore, we need only to show 



the realization of R z (ip) and R y (ip) . To realize the gate 
R y (ip), we let the phases cfro — <fr\ = 0. The U (2) val- 
ued connections can be derived as A^ = —isin(9)<j y} 
and Ag = 0. After an adiabatic evolution along a 
closed path C in the parameter space M, the related 
unitary transformation (holonomy) is just the rotation 
about y axis U (C) = exp [i<p (C) a y ] , where the angle 
ip (C) = — (J c sin (0) d£_) is dependent of the loop C. To 
achieve the gate R z (ip), we can set (fr = £ = (i.e. 
Clo = 0) and change 9 and cfri adiabatically. In this case, 
the non-zero connection is A c f >1 = sin 2 (9) (1 — <j z ). 
As a result, the holonomy associated with a close path 
C can be written as U (C) = e~ ix ^ e ix ^ c ,z where 
X (C) = — \ J sin 2 (9)d<fr\. It is easy to see that U (C) 
is just a rotation about z axis up to a global phase. 

We now illustrate how to realize the controllcd-PHASE 
gate as well as the controlled-NOT gate via the non- 
Abelian holonomy in the present system of SQUID 
qubits, which is the main and novel result of the present 
paper. At this stage, we consider the (3+l)-type energy- 
level structure of two rf-SQUIDs in a microwave cavity. 
We assume that the |e) < — ► \g) transition of each SQUID 
is coupled to the single mode cavity field. Also, we set the 
level spacing oj ei between the states (|e), \at) , i = 0, 1) 
to be different. The transitions |e) < — ► |ao) ( or l a i)) m 
both SQUIDs are coupled to two distinguishable classi- 
cal microwave pulses with different frequencies [15] and 
the " 3-photon resonance" condition of each SQUID is set 
to be satisfied. The Hamiltonian of such SQUIDs in a 
microwave cavity may be written as 



H=J2 [n? ) |e) I <a |+fi? ) |e) J <ai| + ft.c 

1=1,2 



(3) 



+ E [fl (0 « Ie>, <ffl + fe-c] + X) A( ° l e >i < e l ■ 

1=1,2 1=1.2 



Here, A« 



,(') 



,(') 



^ei ~ ^t' 1S tnc dctunning 



,(') 



of the i-th SQUID and Clf^ is the Rabi frequency of the 
microwave pulse coupled to the transition \e) l <-> |oj)j of 
the Z-th SQUID. 

We still choose the states \ai) l i = 0, 1 as the com- 
putational basis of the Z-th qubit. The two SQUIDs are 
coupled indirectly via the single mode cavity field. As in 
the previous discussions for the single-qubit gate case, it 
is seen that when all of the Rabi frequencies Ctf^ of the 
classical microwave pluses are set to zero and the cavity 
is cooled to zero temperature, the two-qubit state, which 
can be written in terms of \cii) 1 \dj) 2 , is isolated from 
the state of the cavity and does not change with time. 
In the following, we show that the holonomy, associated 
with the adiabatic evolution of flf^ along a closed path 
C starting from the point O-^ = (we also refer to this 
point as O) in the parameter space M, can be used to 
achieve a controllcd-PHASE gate or a controlled-NOT 
gate of the two qubits. 
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We first choose Q 



,(i) 



tan [£(')] 



^ and fi« 



ff ( 2 ) tan sec e**i" ( i = 1,2), and take £«, 

9^ and <^ as the independent coordinates in the param- 
eter space M. To realize the controlled-PHASE gate, we 
set #W = 0^ = 0^ = 0. This means that we set the 
Rabi frequencies fij 1 and fi^ to zero and only change 
Oq and f2 Q 2 ^ adiabatically. It is found that the sub- 
space spanned by the states |aj) 1 |oj) 2 |0) c , |a») 1 |<?) 2 |l) c , 
l<7>il<?> 2 |2) c , If?>! I«i> 2 |1> C , |e) 1 |o i ) 2 |0) c , |ai>! |e> 2 |0> c , 
|e)il2> 2 |l> c , lfl>i|e> 2 |l) c and \e) l |e) 2 |0) c = 0, 1) is 
an invariance subspace (we call it as the subspace I) of 
the Hamiltonian (3). After some tedious derivations, the 
dark states, i.e. the basis of the eigenspace of Hamilto- 
nian (3) with zero eigenvalue, can be obtained as: 



Ara ) = Aootsin 


e\ 


sin 






— V2 cos 


> (i) 


sin 




e^ 2) |ao) 1 l5) 2 |l) c 


+V2cos 


> (i) 


cos 


c (2) 


l«o>i |a ) 2 |0) c 


— \/2sm 
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|5)iK) 2 |1) c ] 



|An) - cos |ao>! |oi> 2 |0) c (4) 



sin 



I 

e (2) ] I 

f(2) 
Ol>l l ft l) 2 1°) 



COS 



— sin 



a o)i |ai) 2 |0) c 

|5>iK> 2 |1>c 
«i)i K) 2 |0) c 

^ |oi>i |5) 2 |1) C 



l^io) 



l£ll> 



where A 00 = y2-sin 2 sin 2 [C (2) ] • It is obvious 
that at the point O where £W = £( 2 ' = 0, the above dark 
states are single valued: \Dij (O)) = |ttj) 2 |0) c f° r 
i,j = 0, 1. The non-zero elements of U (4)-valued connec- 
tions arc just A°°f = zA m 2 {2 - sin 2 } sin 2 [^] 

and A 1 ^) = zsin 2 ■ When the system evolves adi- 

abatically along a closed path C in the parameter space 
M and returns to the point O, the associated holonomy 
can be written as U (C) = e <v{c)\a a ) 2 {a a \^(c)\a a a a )(a a a a \ ^ 

Here, the angles r\ (C) and </> (C) are defined as <f) (C) = 
J A - 2 [2 - sin 2 - A 2 ] sin 2 [^] and „ (C) = 

/ sin 2 [^] d^. In the above expression of U (C), 
|a ) 2 (oo| is the projection operator of the second qubit 
to the state |ao) 2 . It is easy to see that the holonomy 
U (C) is the production of a control phase gate operation 
and a single-qubit rotation about the z axis operated on 
the second qubit. When (j> (C) ^ 0, U (C) is an nontrivial 
two-qubit operation. If we choose the path C which sat- 
isfing -q (C) = 0, we can obtain an explicit control phase 
gate via the holonomy U (C) = e i<t>{C)\a a )(a a \ _ 

Moreover, we may also realize the controllcd-NOT by 
setting 6>W = <f>^ = = and choose 9^ 2 \ ^ 



and i^ 2 -* as the control parameters. In this case, the 
dark states can also be obtained with some straight- 
forward derivations, but have quite complicated forms 
(not presented here). Our main result is that, when 
-J c sin [6»( 2 )] d£ (2) -9(C) = tt/2, the holonomy asso- 
ciated with a closed path can be expressed as U (C) = 
e-^4 2) [9 (C)} R { P (f ) Rf ] (f ) UcnRP (-f ), where 



U C n = |a )i (aol ® I (2) + \a 1 } 1 (<ii| <g> ct. 



(2) 



(5) 



is just the controlled-NOT operation. Here, 1^ is the 
identity operator of the sec- 

ond qubit and the angle 9 (C) are defined as 9 (C) = 

- J c /T 1 A- 1 1 A 00 3 A 2 ^ (2) - Sc P' 1 A oi lA oo 3A i^ (1) . where 
A o is defined as before and the other coefficients are 

defined as Aqi = 



(i)l 



(2)1 



sm' 



= - (1/2) A^ 1 A^ 1 sin 2 sin [2^] sin [9^ 



[f< 2 >], 
'1 and 



/3 = \/l + a 2 . Therefore, up to a global phase factor, 
U (C) is just the product of control not operation and 
some single-qubit rotations. 

In most cases, the Hamiltonian (3) has a four- 
dimension engenspace with zero eigenvalue in the invari- 
ant subspace I. The basis of the eigenspace are just the 
dark states |-DjA. Nevertheless, it is pointed out that, in 
some very special cases, there may be accidental degen- 
eracy in a sub-manifold (we call it the AD sub-manifold) 
of the parameter space M. In the AD sub-manifold, in 
addition to the dark states \Dij), the Hamiltonian has 
another two eigenstates with zero eigenvalue and thus 
the dimension of the Hamiltonian's eigenspace with zero 
eigenvalue is six rather than four. It is apparently that 
if the path of the adiabatic evolution of the Rabi fre- 
quencies cross the AD manifold, there might be a tran- 
sition from the four dark states \Dij) to the two exter- 
nal states. To avoid this kind of unwanted transition, 
we should control the evolution path of the Rabi fre- 
quencies in the parameter space M be far away enough 
from the AD sub-manifold. On the other hand, since 
the evolution of the Rabi frequencies in M is assumed 
to begin and end at the same point O where all the 
Rabi frequencies are set to zero, the accidental degen- 
eracy at O should be avoided. This can be implemented 
by adjusting the coupling strength g^> via controlling 
the position of the SQUIDs in the cavity, or the detun- 
ing AW of each SQUID via changing the bias fields. For 
instance, if the conditions = A( 2 ' ^ 0,g^ = or 
A(i) _ ^(2) = o, 5 ( 1 ) ^ g( 2 ) are satisfied, the accidental 
degeneracy at point O can be avoided. 

Since the dark states have a non-zero projection to 
the single or two photon states, the photon dissipation 
caused by the imperfection would tend to destroy the 
dark states. To evaluate the influence of the photon dis- 
sipation, the Schorcdinger equation controlled by the ef- 
fective Hamiltonian 



H eff =H(t)-i {n/2)a}a 



(6) 
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needs to be solved. Here, H (t) is the Hamiltonian defined 
in Eq. (3) (or Eq. (2)), and k the photon decay rate. 
The dissipation term —i (k/2) a^a can be considered as a 
perturbation. As a result of the first order perturbation 
theory, the fidelity of the quantum gate operation, i.e., 
the probability of the ideal finial state, may be expressed 
as 

F « 1 - k f (n) ph dt (7) 

JQ 

where T is the operation time and (n) ph — 
(ty (t) \ a^a |\I> (t)) is the instantaneous expectation value 
of photon number. Therefore, the condition under which 
the influence of photon dissipation can be neglected is 
simply k J Q (n) h dt << 1. On the other hand, since 
our scheme is based on the adiabatic evolution of the 
quantum states, the adiabatic condition should be sat- 
isfied, which can be expressed as f2T >> 1, where £7 is 
the energy gap between the dark state and other eigen- 
states of the Hamiltonian [5]. Here, f2 has the same order 
of amplitude as the SQUID-cavity coupling constant g. 
In practical quantum gate operations, we always have 
T ~ 10 3 ff -1 . Then the condition k (n) ph dt « 1 
can be satisfied when g/n > 10 4 . The coupling con- 
stant of the SQUID and the cavity available at present 
is g ~ 1.8 x lifs- 1 [15]. The high quality factor of the 
cavity Q — 10 6 ~ 10 s might be achieved experimentally 
[22]. This will lead to k J q T (n) ph dt < 1CT 2 and thus the 
fidelity F ~ 1. 




r. 



FIG. 2. The fidelities of a control phase gate for the ini- 
tial states \ao} 1 |oo) 2 (solid line), \a ) 1 |oi) 2 (dashed line), 
\ai) 1 |ao) 2 (dotted line), and |ai) 1 |ai) 2 (dashed-dotted line), 
where the quantity n is defined by the relation g/n = 10". 

Finally, let us look into in some detail a typical con- 
trol phase gate operation discussed before. In this oper- 
ation, we assume g^ — 2g^ = g = 1.8 x 10 8 s _1 and 
A^ 1 ) = A^ 2 ) = 0. The amplitudes of the Rabi frequencies 
and ^ are varied following the Gaussian functions 
of time: n { 1} = 2.5ge-(^) 2 ,^ 2) = ge'^^)' ' e^o 
where r = 144g _1 . The phase 4>^ (t) is set to be 
a hyperbolic tangent function of time: (j) (t) = n x 
[l + tanh ( Q y 5r )] ■ As in the above discussions, the con- 
trol gate operation can be written as an unitary trans- 



formation U = e *'?l a o> 2 <oo| e i^|oooo>(oooo| ) w here 4> w tt/6 
and r\ w 4. The operation time of this gate is about 
8 x 10 2 <7 _1 . We estimate the fidelity of this operation 
using Eq. (7). In Fig. 2, the fidelity of the quantum 
gate operation with four possible initial states are plot- 
ted as a function of the ratio g/n. It is seen that when 
gj ' k <~ 10 3 — 10 4 , the fidelity is larger than 90%. In par- 
ticular, whenever g/n <~ 10 5 , the fidelity is improved to 
reach 99%. 
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